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Abstract 

We have developed the full theory of a synchronously pumped type I optical parametric oscillator (SPOPO). We derive 
expressions for the oscillation threshold and the characteristics of the generated mode-locked signal beam. We calculate the 
output quantum fluctuations of the device, and find that, in the degenerate case (coincident signal and idler set of frequencies), 
perfect squeezing is obtained when one approaches threshold from below for a well defined "super-mode', or frequency comb, 
consisting of a coherent linear superposition of signal modes of different frequencies which are resonant in the cavity. 

PACS numbers: 42.50.Dv, 42.65. Yj, 42.65.Re 



Optical Parametric Oscillators are among the best 
sources of squeezed 0, correlated and entangled Q 
light in the so-called continuous variable regime. They 
have allowed physicists to successfully implement demon- 
stration experiments for high sensitivity optical measure- 
ments and quantum information protocols. In order to 
maximize the quantum effects, one needs to optimize 
the parametric down-conversion process. This has been 
achieved so far by using either intense pump lasers or res- 
onant cavities. Having in mind that the parametric pro- 
cess is an almost instantaneous one, femtosecond mode- 
locked lasers are the best pump sources in this respect, 
as they generate very high peak optical powers with high 
coherence properties. Furthermore, they minimize the 
thermal effects in the linear crystal which often hamper 
the normal operation of parametric devices. Mode-locked 
lasers have been already used extensively to generate non 
classical light, either to pump a parametric crystal or 
an optical fiber However in such single-path config- 
urations, perfect quantum properties are only obtained 
when the pump power goes to infinity. This is the reason 
why mode-locking is often associated to Q-switching and 
pulse amplification 6] in order to reach even higher peak 
powers, at the expense of a loss in the coherence prop- 
erties between the successive pump pulses. In contrast, 
intracavity devices produce perfect quantum properties 
for a finite power, namely the oscillation threshold of the 
device. It is therefore tempting to consider devices in 
which one takes advantage of the beneficial effects of both 
high peak powers and resonant cavity build-up. Such de- 
vices exist: they are the so-called synchronously pumped 
OPOs or SPOPOs. In a SPOPO, the cavity round-trip 
time is equal to the repetition rate of the mode-locked 
laser, so that the effect of the successive intense pump 
pulses add coherently, thus reducing considerably its os- 
cillation threshold. Such SPOPOs have already been 
implemented as efficient sources of tunable ultra-short 
pulses 0, @, El EJ El E3 an d their temporal proper- 
ties have been theoretically investigated [lifl, LLJ, ■ To 
the best of our knowledge, they have never been used 
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FIG. 1: Synchronously pumped OPO 



so far to generate genuine quantum effects. Perform- 
ing a quantum analysis of this device, we theoretically 
show in this paper that perfect squeezing can indeed 
be obtained in SPOPOs as in cw OPOs. The squeezed 
mode is not a usual single frequency mode, but instead 
a "super-mode", which is a well defined linear combi- 
nation of signal modes of different frequencies, forming 
the frequency comb that will oscillate above the SPOPO 
oscillation threshold. 

Let us first precise the model that we use (figure [Q. 
We consider a ring cavity of optical length L containing 
a type I parametric crystal of thickness I, Degenerate 
phase matching is assumed, what means that the phase- 
matching condition is fulfilled for frequencies 2wo and too. 
This amounts to saying that n (2loq) = n (u>o) = no, n (u>) 
being the refractive index of the crystal at frequency ui. 
The mode-locked pump laser, having a repetition rate 
n/2ir = c/L, is tuned so that the frequency of one of its 
modes is equal to 2lo . The electric field generated by 
the pump mode-locked laser can be expressed as: 
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where P is the average laser power per unit area, a m 
the normalized (%2 m \ a m\ 2 = 1) complex spectral com- 
ponent of longitudinal mode labelled by the integer index 
to, and m — corresponds to the phase-matched mode. 
For the sake of simplicity in this first approach to the 
problem, we will take the modal coefficients a m as real 
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numbers, thus excluding chirped pump pulses. As al- 
ready mentioned, the SPOPO cavity length is adjusted 
so that its free spectral range coincides with that of the 
pumping laser. In the nonlinear crystal, pump photons 
belonging to all the different longitudinal pump modes 
are converted into signal and idler photons via the para- 
metric interaction. In addition we will assume here that 
we are in the ideal case of doubly resonant degenerate 
operation, meaning that the among all the OPO cavity 
resonant frequencies, they are all the pump mode fre- 
quencies uj p ^ m = 2luq + m£l but also all the frequencies 
u> Sl q = ujq + qfl around the phase-matched subharmonic 
frequency luq. The intracavity electric field generated by 
the parametric interaction will then be a superposition 
of fields oscillating at frequencies oj s>q . We will finally 
call 7 P and 7 S , the cavity damping rates for the pump 
and signal modes. Note that the free spectral range f2 is 
assumed to be the same in the pump and in the signal 
spectral regions. This is necessary for an efficient in- 
tracavity parametric down conversion and requires, from 
the experimental viewpoint, the use of extra dispersive 
elements inside the cavity that compensate for the dis- 
persion of the crystal. At the quantum level, the signal 
field, taken at the middle of the crystal, is represented 
by the quantum operator E s which can be written as: 



E s (t) = } j i£ s , q s q {t)e 
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where s q are the annihilation operators for the q signal 
mode in the interaction picture. £ s _ q is the single photon 
field amplitude, equal to yJhw SA j2e^n (uj s , q ) AL, and A 
its effective transverse area. The following Heisenberg 
equations for the field operators can be derived using the 
standard methods. The detail of the derivation will be 
given in a forthcoming publication. Below threshold, and 
in the linearized regime for the pump fluctuations, they 
read: 

ds m 

~dT 

where a is the normalized pump amplitude 
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in which Pq is the single mode c.w. oscillation threshold: 
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with x the crystal nonlinear susceptibility. £ m , q is the 
product of a phase-mismatch factor by the pump spectral 
normalized amplitude a rn + q : 



Iy m,q 



The phase mismatch angle 
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can be computed using a Taylor expansion around 2wo 
for the pump wave vectors k p ^ m and around ujq for the 
signal wave vectors k s q : 



wq ~ f3i (to + q) + (3 2p (m + q f - (3 2s (m 2 + q 2 ) , 

re ft = in (k' p - k' s ) i, p 2p = \n 2 k;i, (3 2s = \n 2 
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where 1 = Afi (k' p - k' s ) I, (3 2p = \n 2 k' p % f3 2s = \n 2 k'^l. 
k! and k" are the first and second derivative of the wave 
vector with respect to frequency. Finally Si n ,m are the 
input signal field operators at frequency u) s ,m transmit- 
ted through the coupling mirror. When the input is the 
vacuum state, which we consider here, their only non-null 
correlations are: 
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In order to get Eqs. ©, we had to assume that £ s ^ m ~ 
£ Si q for all m and to neglect the dispersion of the non- 
linear susceptibility. Both approximations require that 
pulses are not too short. In usual practical conditions, 
pulse durations should be longer than 20-30 fs. Let us 
first determine the average values of the generated fields. 
They are determined by the "classical" counterpart of 
Eq. ©, removing the input noise terms, and replacing 
the operators by complex numbers. The solution of these 



equations is of the form s m (t) = Sk 



where k is an 



index labelling the different solutions. The parameters 
Sk m an d Afc obey the following eigenvalue equation: 
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As matrix C is both self-adjoint and real (C m , q — £ q .m 
real, see Eqs. its eigenvalues A fc and eigenvec- 

tors -Lfc, of components L^ m , are all real. As 7„ and cr are 
also real, there exist two sets of solutions of Eqs. QlOp. 

The first set is given by 



that we will call and S { k 2 
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corresponding eigenvalues: 
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Let us now label by index k = the solution of maximum 
value of |Afc|: |Ao| = max{|Afe|}. When a |Ao| is smaller 
than 1, all the rates A^ are negative, which implies that 
the null solution for the steady state signal field is stable. 
For the simplicity of notations, we will take Ao positive 
in the following 16] . The SPOPO reaches its oscillation 
threshold when a takes the value 1/Aq, i.e. for a pump 
power P = P t hr equal to: 



P thr = Po/Al 



(12) 



We can now define the normalized amplitude pumping 
rate r by r — a Ao, so that threshold occurs at r = 1. 
We will call eigen-spectrum the set of Sk, m values for 
a given k, which corresponds physically to the differ- 
ent spectral components of the signal field, and critical 
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eigen-spectrum Sp*^, the one associated with Aq 1 ~ /I , which 
changes sign at threshold. Above threshold, this criti- 
cal mode will be the "lasing" one, i.e. the one having 
a non-zero mean amplitude when r > 1. Let us note 
that the eigen-spectrum in quadrature with respect to 
the critical one, Sq^ = iSq , has an associated eigen- 
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value A 



(-) _ 



-27,5 at threshold. Furthermore, equation 



implies that all the damping rates a! are com- 
prised below threshold between — 2-f s and 0, and that, 
whatever the pump intensity, all the eigenvalues At (r) 
lie between Aq (r) and Aq ' (r). These properties will 
be useful for the study of squeezing. To determine the 
threshold, we must find Ao, and therefore diagonalizc 
C, which cannot be done analytically in the general case. 
The detailed study of the SPOPO threshold will be made 
in a forthcoming publication. Here we will consider a 
special case that leads to simple calculations and corre- 
sponds to an optimized situation. We assume first that 
0x = (equal group velocities) and /?2 S = 2/?2 P (matched 
group velocity dispersions). We take simple forms for the 
phase matching coefficient 



Cr. 
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and for the normalized pump spectrum: 
a m = 7 r- 1 / 4 A; 1 /2 e -^ 
If in addition we impose the relation 
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giving comparable widths for the phase matching curve 
and the laser spectrum, we can analytically diagonalize 
the matrix C rn . q . We find that its (unnormalized) eigen- 
vectors have the following components: 
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Hk being the Hermite polynomial of order k, and that 
its eigenvalues are: 



is very fast, the OPO actually oscillates when the in- 
stantaneous peak power of the pump exceeds the c.w. 
threshold. Having in mind that Pq is on the order of a 
few tens of mW and that the number of coherently os- 
cillating modes in a mode-locked laser can easily exceed 
10 4 , one can then envision SPOPO thresholds of the or- 
der of a microWatt. Experimental implementations of 
such SPOPOs have not reached such ultra-low thresh- 
olds because of the difficulty to fulfill simultaneously, in 
a real experiment, all these precise conditions. We can 
now determine the squeezing properties of the signal field 
in a SPOPO below threshold. This is done by using the 
SPOPO linearized quantum equations. Let us introduce 
the operator Si n ^{£) by: 



E. 
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Sin,k(t),S! n>k (t') 



Because J2 m \ L k,m\ 2 = 1, one has 

S(t — t'): §i n f. is the annihilation operator of a com- 
bination of modes of different frequencies, which are the 
eigen-modes of the linearized evolution equation J3J). The 
corresponding creation operator applied to vacuum state 
creates a photon in a single mode, which can be labelled 
as " super-mode" , which globally describes a frequency 
comb. Defining in the same way as in (|19|l the intracav- 
ity operator Sk(t), one can then write: 



dt 



-lsS k + -fsvAkSl + ^/WsS^k, (20) 



Let us now define quadrature hcrmitian operators Sj. 
by: 



(±) 



(16) which obey the following equations: 
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(22) 
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A fc = 5 
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This gives finally a SPOPO threshold (Eq. {Hty) equal 
to: 
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2P 
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the associated eigenspectrum being a Gaussian of 
width A p . Relation 118fl means that the usual c.w. OPO 
oscillation threshold is divided by a quantity roughly 
equal to the number of pump modes, at least in the fully 
optimized SPOPO that we have considered here. This 
has a simple explanation: since the parametric coupling 



with At gi ven by Eq. 111|) . These relations enable us 
to determine the intracavity quadrature operators in the 
Fourier domain (to) 



(±>/ 



(24) 



Finally, using the usual input-output relation on the 
coupling mirror: 



(25) 



which extends by linearity to any super-mode operator 
because the mirror is assumed to have a transmission 
independent of the mode frequency, one obtains finally 
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the following expression for the output signal super-mode 
quadrature component in Fourier space: 

5 ° ut < fc(a;) " ^T±7AVA^T^ 5i ^ fc(a;) ' (26) 

These expressions are particularly simple for the critical 
mode quadrature components (k = 0): 



7,(l±r)-tw x+ 



-7s (1 =p r) + 
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The variance of these operators can be indeed mea- 
sured using the usual balanced homodyne detection 
scheme: the local oscillator is in the present case a coher- 
ent mode-locked multimode field E^, (t) having the same 
repetition rate as the pump laser: 



(28) 



where J2m \ e m\ 2 = 1, and cl is the local oscillator field 
total amplitude factor. Assuming that the photodetec- 
tors measure the intensity of the Fourier components of 
the photocurrent averaged over many successive pulses, 
the balanced homodyne detection scheme measures the 
variance of the fluctuations of the projection of the out- 
put field on the local oscillator mode when the mean field 
generated by the OPO is zero, which is the case below 
threshold. As a result, when the coefficients e m of the 
local oscillator field spectral decomposition are equal to 
the coefficients Lk t m of the fc-th super-mode, one mea- 
sures the two following variances, depending on the local 
oscillator phase: 

Vn^ 2(1 ' rAfc/A0) ^ 2 (29) 
k 7s 2 (l + rA fc /A ) 2 +c 2 

V+ (u;) = V k - (uj)- 1 (30) 



Eqs. 



shows that quantum noise reduction below 



the standard quantum limit (equal here to 1) is achieved 
for any super-mode characterized by a non-zero A& value 
and that the smallest fluctuations are obtained close to 
threshold and at zero Fourier frequency: 



A - [Afcj 
A + |A fc | 



(31) 



In particular, if one uses as the local oscillator the critical 
mode k = 0, identical to the one oscillating just above 
the threshold r = 1, one then gets perfect squeezing just 
below threshold and at zero noise frequency, just like in 
the c.w. single mode case. But modes of k ^ may be 
also significantly squeezed, provided that |Afc/Ao| is not 
much different from 1. In conclusion, we have studied 
the quantum behaviour of a degenerate synchronously- 
pumped OPO, which seems at first sight a highly 
multimode system, since it involves roughly 10 5 different 



usual single frequency modes for a 100/s pulse. We have 
shown that its properties are more easily understood if 
one considers the "super-modes", linear combinations 
of all these modes that are eigen-modes of the SPOPO 
set of evolution equations and describe in a global way 
the frequency comb -or, equivalently, the train of pulses- 
generated by the SPOPO. The super-mode of minimum 
threshold plays a particular role, as it is the one which 
turns out to be perfectly squeezed at threshold and will 
oscillate above threshold. The present paper gives a first 
example of the high interest of studying frequency combs 
at the quantum level, as they merge the advantages of 
two already well-known non-classical states of light: the 
c.w. light beams, with their high degree of coherence 
and reproducibility and the single pulses of light, with 
their high peak power enhancing the non-linear effects 
necessary to produce pure quantum effects. 
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